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Feedback Linearization Continued 


Recall “strict feedback systems" discussed in Lecture 13: 
ty = fil%r) + 91 (41) x2 


Xo = fa(x1, x2) + 92(x1,%2)x3 


X3 = fa(x1,X2,%3) + g3(x1,%2,%3) x4 (1) 


Xn = fn(x) + Sn(x)u. 
Such systems are feedback linearizable when ¢;(x1,...,x;) # 0 near 
the origin, i = 1,2,--- ,n, because the relative degree is n with the 
choice of output y = h(x) = x1: 
y) = Lth(x) + gr (x1) g2 (x12) --- Bn(x) 
: eS Ey EE 
LoL? 'h(x) #0 


Feedback linearizability is lost when g;(0) = 0 for some i; however, 
backstepping may be applicable as the following example illustrates: 


Example 1: 
a x t0 
Xo = U. 
Treat x2 as virtual control and let #1(x,) = —x, which stabilizes the 


x1-subsystem, as verified with Lyapunov function Vj (x1) = 4x2. 


Then Z2 := X2 — &1(x}) satisfies Z7 = u — &1, and 


Vi 
—— xt — kyz0 — x2X9 x kp (x2 + x1) 


cm ae 


achieves global asymptotic stability: 


Za V = —xy4 — kgzs. 


In contrast the system is not feedback linearizable, because condition 
(C1) in the theorem for feedback linearizability (Lecture 20, p.4) fails. 
To see this note that 


2 


pee) = |, w= [1], agate) =r) = [9] 
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thus, with n = 2 and xp = 0, 


[g(x0) adyg(t0) «.. ad? *g(x0)] = f | 


which is rank deficient. 


Multi-Input Multi-Output Systems 


Consider now a MIMO system with m inputs and m outputs: 


m 


e = Fa) eae (2) 
i=l 
yi = A(x), i=1,---,m. 


Let r; denote the number of times we need to differentiate y; to hit at 
least one input. Then, 


yO] PLB (x) ] [bg kM (x) e+ Lgl M(x) ] fay 

ol ee les 

vn Lehn (x) Lg, Le him (x) — Len, Li him (x) Um 
=: B(x) =a A(x) 


If A(x) is nonsingular, then the feedback law 
u = A(x)"(-B(x) +0) 


input/output linearizes the system, creating m decoupled chains of 
integrators: 
yi) =O, t= Lpes.jm: 


We say that the system has vector relative degree {r1,--+ ,tm} if the 
matrix A(x) defined above is nonsingular. 


Example 2: The kinematic model of a unicycle, depicted below, is 


[4% cos X3 0 
xo| = |sinx3| uy + |0]} Uo, 
Is] Po fy 


where wu is the speed and up is the angular velocity. 
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Let y, = x1 and y2 = x2, and note that 


Yi| — |cosx3 Of} |uy 
y2 ~ | sin x3 Of} |u2 ; 
a 


=: A(x) 


Since A(x) is singular, the system does not have a well-defined vector 


relative degree. 
The notion of zero dynamics and the normal form can be extended to 

MIMO systems?. If the system has vector relative degree {r1,--- ,1m}, * see, e.¢., Sastry, Section 9.3 
then r:= 77 +--+ +%m <n and 


G2 [fa (x) Lyla (x) ++ LET (x) +++ B(x) Lyle (x) ++ Lm (x)] 


defines a partial set of coordinates. As in normal form discussed in 
Lecture 19, one can find n — r additional functions z1(x),--+ ,Zn—r(X) 
so that x ++ (z,C) is a complete coordinate transformation. 


Full-state feedback linearization amounts to finding m output func- 

tions /,--- , Hm such that the system has vector relative degree 

{11,--+ tm} with 7; +---+1%m = n. Necessary and sufficient con- 

ditions for the existence of such functions, analogous to those in 

Lecture 20 for SISO systems, are available3. 3 see, e.g., Sastry, Proposition 9.16 


Example 3: Consider the following model of a planar vertical take-off 


and landing (PVTOL) aircraft* 4 Sastry, Section 10.4.2 
¥ = —sin(@)u, + pcos(@)u2 
Z = cos(@)u, +psin(@)u2 —1 
6 = up, 


where p/ is a constant that accounts for the coupling between the 
rolling moment and translational acceleration, and —1 in the second 
equation is the gravitational acceleration, normalized to unity by 
appropriately scaling the variables. 


N 
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If we take x and z as the two outputs we get 


x) | 0 4 —sin@ pcos} | uy 
Zz —1 cos@ =psin®@| |u2 
A(8) 
where A(@) is invertible when p # 0: 


=cos@ -=siné 


—sin@ cosé@ 
# # 


A (G)= F 1 


Thus the systems has vector relative degree {2,2} when p # 0, and 
the input/output linearizing controller is 


uy —sin@ cosé 0 V1 
=}1 0 le Q + . 
uy i COS 1 Sin 1 V2 


The zero dynamics is obtained by substituting u; = i sin 0, needed 
to maintain z at a constant value and Z at zero, in the dynamical 
equation for 0: 


es 1 
6 = —siné. 
H 


The system is nonminimum phase for pt > 0, since 6 = 0 is unstable. 


Drift-Free Systems 


Suppose f(x) = 0 for all x in (2). Such system are called drift-free and 
encompass linear systems of the form 


x= Bu, x €R", uc R". 


Assuming the columns of the n x m matrix B are linearly indepen- 
dent, we can find n — m row vectors T;,i = 1,--- ,m—m, such that 


T;B =0. 


This means that #;(x) := T;x satisfies 


THlx()) =0 > gi(x(t)) = gu(o(0)) ) 


regardless of the control inputs. Since there are n — m such con- 
straints, controllability is not possible in drift-free linear systems with 
fewer control inputs than the state dimension (m < n). 


The Frobenius Theorem (Lecture 20) implies that constraints of the 
form (3), called holonomic constraints, also exist for nonlinear drift-free 
systems 


x= gil) ui (4) 


EE C222/ME C237 - SPRING’18 - LECTURE 21 NOTES 


when the distribution A = span{g1,--- ,2m} is nonsingular and 
involutive. 


When A is non-involutive, however, controllability may be possible 
with m <n — another essentially nonlinear phenomenon. 


Indeed, Chow’s Theorem states that (4) is controllable if the involutive 

closure> of A = span{g1,--- , 8m} has dimension n. This condition 5 the smallest involutive distribution 
means that the Lie brackets of g1,--- , %m Span new dimensions that that containts A 
are not already spanned by these basis vector fields. Drift-free sys- 


tems satisfying Chow’s Theorem are called nonholonomic. 


Example 4: Recall the unicycle model discussed in Example 2, where 


COS X3 0 — sin x3 
gi(x) = |sinx3], g2= 0], and [gi,go](x) = | cosxs 
0 1 0 


A = span{g1,82} is non-involutive, as [g1,%2] generates a new di- 
rection. Taken together, 91, g2, and [21,22] span the entire three- 
dimensional space at each point x; therefore, the system is control- 
lable by Chow’s Theorem. This conclusion sheds light on how paral- 


lel parking is possible despite lack of sideways actuation. 


To present an interpretation of the Lie bracket [¢1,g2], we let ®¥' (xq) 
denote the solution of the system x = g;(x) at time t from initial 
condition x9. Then it can be shown that 


d, ?(, 8 (GP? (GF (x0)))) = PF lg1-gal(x0) + O(F), 


which suggests that motion in the direction of the Lie bracket [¢1, 99] 
can be generated by alternating actuation of the two inputs u, and u 
with positive and negative signs, as one does in parallel parking. 


